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We examine the validity of the recently proposed semi- 
Poisson level spacing distribution function P{S), which char- 
acterizes 'critical quantum chaos', in 2D disordered systems 
with spin-orbit coupling. At the Anderson transition we show 
that the semi-Poisson P{S) can describe closely the criti- 
cal distribution obtained with averaged boundary conditions, 
over Dirichlet in one direction with periodic in the other and 
Dirichlet in both directions. We also obtain a sub-Poisson 
linear number variance T,2{E) ~ xo + X^i with asymptotic 
value X ~ 0.07. The obtained critical statistics, intermedi- 
ate between Wigner and Poisson, is relevant for disordered 
systems and chaotic models. 



In mesoscopic physics the effect of disorder on the elec- 
tron propagation leads to the zero-temperature quantum 
Anderson metal-insulator transition, which arises from 
the competition between quantum tuncUing and inter- 
ference, as a function of disorder. For weak disorder the 
electrons diffuse, due to quantum tunelling, and the sys- 
tem is metallic with correlated chaotic energy levels and 
'level-repulsion' described by Wigner statistics In 
the case of strong disorder the electrons localize in ran- 
dom positions, due to quantum interference, and the sys- 
tem becomes insulating, having non-chaotic completely 
uncorrelated (random) energy levels which show 'level at- 
traction' and obey ordinary Poisson statistics. In order 
to see the metal-insulator transition high enough space 
dimensionality (usually greater than 2) is required and 
at the critical point, which corresponds to an interme- 
diate value of disorder, the level statistics changes from 
Wigner to Poisson The critical electrons are nei- 

ther extended nor localized and it is believed that a new 
universal critical statistics, intermediate between Wigner 
and Poisson, should apply [§|J|]. We aim to address the 
question of the critical statictics in two dimensions (2-D), 
where a metal-insulator transition occurs in the presence 
of spin-orbit coupling |6). 

The stationary energy levels of 2D quantum billiards 
(e.g. in the form of the stadium), with zero potential 
inside and infinity outside, can also display quantum 
chaotic behavior |Q. The analogies in the level statis- 
tical description bring together the two fields of meso- 
scopic physics and quantum chaos and have been ex- 
ploited in the past for understanding important phenom- 
ena in both areas. In this respect, Wigner statistics was 
originally conjectured to apply for quantum systems with 
chaotic classical dynamics, since the levels resemble the 



eigenvalues found in appropriate random matrix ensem- 
bles, introduced long ago 0,^. On the other hand, in- 
tegrable systems correspond to Poisson statistics having 
completely uncorrelated (random) eigenvalues. The key 
question is again what happens at criticality, between 
chaos and integrability, similarly to the transition be- 
tween metal and insulator. Recently, a new distribu- 
tion was proposed to describe critical levels statistics 
[^,|o|, which contains both Wigner and Poisson features, 
as the main theme of what is called 'critical quantum 
chaos'. This intermediate distribution can be derived 
from a short range plasma model j|] and was realized 
in pseudo-integrable systems, such as the classically non- 
integrable but of zero metric entropy rational triangle 
billiards , and corresponds to other solvable models 
[|l2| , ^ . In disordered systems the intermediate distri- 
bution, named semi-Poisson, was shown to characterize 
critical states at the iD metal-insulator transition |l^] 
and the energy levels of few electrons in the presence of 
interactions Itl4|. 
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FIG. 1. Energy levels E versus disorder W , for a 2D disordered 
system with spin-orbit coupling. The metal-insulator critical point 
Wc = 8.55, marked with a broken line, separates chaotic levels 
(on the left) with 'level-repulsion' and 'spectral rigidity' described 
by Wigner statistics, from non-chaotic levels (on the right) with 
'level attraction' and 'spectral randomness' described by Poisson 
statistics. In this paper we address what happens at the metal- 
insulator transition (broken line). 

In Fig.l the energy levels E obtained from our 2D dis- 



1 



ordered system are displayed as a function of disorder W. 
The behavior of the levels is seen clearly to change at the 
critical point of the metal-insulator transition Wc, which 
separates chaotic levels (on the left) from non-chaotic 
levels (on the right). The chaotic levels are more reg- 
ular (correlated) than the non-chaotic levels, which arc 
uncorrelated (random). The 'level- repulsion' effect can 
be seen on the chaotic levels and the 'level-attraction' 
on the non-chaotic levels, where degeneracies exist. The 
displayed levels in Fig. 1 are obtained for one random 
configuration from a system of linear size L = 30, by 
keeping all the energies at the middle of the spectrum, 
making the density of states constant at every W. The 
overall 'spectral rigidity', seen in the metallic chaotic lev- 
els, can be contrasted with the 'spectral randomness' of 
the insulating non-chaotic levels. The surprising result, 
also visualized in this figure, is that the chaotic levels 
appear more 'regular' than the non-chaotic ones. 

The question addressed in this paper is: 'what is the 
level statistics at the critical broken line of Fig. 1?' 
This study is done in connection to the scenario of 'crit- 
ical quantum chaos', which is summarized in a semi- 
Poisson level spacing distribution function P{S) and a 
sub-Poisson linear number variance T,2{E), which mea- 
sures level-fluctuations in a given energy window E. 
Moreover, in order to examine the validity of the semi- 
Poisson P{S) different boundary conditions (BC) must 
be considered, in the spirit of recent important findings 
I pof . The fact that the critical level fluctuations arc at 
the same time scale-invariant (size-independent) and de- 
pendent on BC was explained by invoking the concept of 
the critical conductance . The influence of boundary 
conditions on the critical level statistics, without affect- 
ing the critical disorder, has also been demonstrated for 
a different critical 2D model jl^. 

The main features of 'critical quantum chaos' for sys- 
tems characterized by the universality index P = 1,2,4 
are summarized in: (i) The semi-Poisson P{S) level spac- 
ing distribution which shows Wigner-like repulsion ~ 
at small spacings S << 1 and is exponential, Poisson- 
like, ~ exp(— (/9-I- 1)5*) at large spacings S » 1, overall 
described by the scale-invariant normalized semi-Poisson 
curve 

P{S) = AS^exp{-i(3+l)S), (1) 

with the constant values A ^ A, 27/2 and 3125/24 ob- 
tained from normalization, respectively. The spacing dis- 
tribution P{S), by applying a 'level-unfolding' procedure 
which keeps the level-density constant, corresponds to 
(S) = 1. (ii) The sub-Poisson number variance, which 
defines the level number fluctuations in an energy win- 
dow E, with the mean number proportional to E after 
'unfolding', according to this scenario is 

^2iE)^Xo + xE, (2) 



defining the level compressibility x- The value of x ranges 
between (chaos) and 1 (integrability) and was related 
to the multifractality of the critical wavefunctions [|l6|. 

The considered disordered system displays a transition 
in 2D with energy levels which obey Wigner statistics for 
the metal (with /3 = 4) and Poisson statistics for the insu- 
lator (see Fig. 1). At criticality, where one expects 'crit- 
ical quantum chaos' to apply, numerical work suggested 
level- repulsion in 3D for small S [Q,^ , also later shown in 
2D In order to study carefully the level fluctuations 
in the critical region it is important to identify the cru- 
cial role of BC |jl^ . We find that for the three considered 
kinds of BC the critical distribution function shows level 
repulsion at small spacings and is Poisson-like at large 
spacings. However, when considering an averaged distri- 
bution over the cases: 1) Dirichlet BC in both directions 
and 2) periodic BC in one direction and Dirichlet in the 
other, the obtained distribution is seen to be remarkably 
close to the scale-invariant semi-Poisson curve of Eq. (1), 
appropriate for /? = 4 (see Fig. 3 below). 

The theoretical framework to study the Anderson tran- 
sition can classify disordered systems into three uni- 
versality classes, depending on whether the Hamilto- 
nian preserves the time-reversal invariance or the rota- 
tional invariance, in direct analogy with the random ma- 
trix theory description of quantum chaotic systems . 
Zero spin-orbit corresponds to the orthogonal universal- 
ity class {(3—1) and finite spin-orbit breaks the rota- 
tional invariance so that one obtains the symplectic uni- 
versality class (/3 = 4). In our calculations we consider a 
two-dimensional disordered system, with spin-orbit cou- 
pling for spin-i particles, described by the Hamiltonian 

i 

n = 5]e.c+,c,,, + <^^j,y ' (3) 

where i labels the square lattice sites and a = ±1/2 
is the spin index on each site. The second sum is taken 
over all nearest neighbour lattice pairs and the ran- 
dom on site potential is a spin independent uniformly 
distributed random variable, chosen from a probability 
distribution of width W. In this case the nearest neigh- 
bour hoppings Vij are random 2x2 matrices describing 
spin rotation, due to spin-orbit, on every lattice bond 
(i, j). In the spinor space they are represented by 

-^yy + ifiV'-" 1 - ifiv J ^. ' ^ ' 

where /i denotes the spin-orbit coupling and the V^, 
and V^, defined for every bond are real and inde- 

pendent random variables chosen from a uniform proba- 
bility distribution on the interval [— ^,-1-^]. For the rest 
the spin-orbit strength is fixed to /i = 2 and the disorder 
is chosen to lie exactly at the critical point Wc = 8.55 
§■ 
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We compute the eigenvalues from Eq. (3) by diagonal- 
izing numerically the corresponding Hamiltonian matri- 
ces for large square lattices. The statistical analysis of 
energy levels must be done on a constant density of states 
using an 'unfolding procedure'. In order to achieve the 
level unfolding for the considered disordered system it 
is sufficient to obtain the average of the integrated den- 
sity of states A/", locally at E, by repeating many times 
the disorder configuration creating a statistical ensemble. 
Then the 'raw' spacings Ai ^ Ei ~ -E^-i are replaced by 
the 'unfolded' new ones Si = Mav{Ei) — Nav{Ei-i) « 
{E, - E,_i f^l^-^^^ = {E, - E,_i)/A, where A is the 
local mean spacing around Ei or equivalently the inverse 
density of states of the raw data. In the numerical calcu- 
lations we considered eigenvalues within the energy win- 
dow [-2, 2] performing 2700, 1200, 675 and 794 random 
configuration runs in each case, with L = 20, 30, 40 and 
60, respectively. The total number of eigenvalues from 
all random configurations for each BC is about 400,000 
for L = 20, 30, 40 and 1, 000, 000 for L = 60. These 'raw' 
data were 'unfolded' in the described way. 




FIG. 2. This figure shows the variation of the critical P{S) 
for three choices of BC for a system of hnear size L. The mean 
distribution over the three cases is shown by the continuous black 
line. The Wigner distribution (dotted line) is also shown. 

The obtained P{S) at criticality is shown in Fig. 2 for 
the three different choices of BC. The computed curves 
are, clearly, very different, in agreement with the corre- 
sponding 3D results for P = 1 They are very dif- 
ferent from the Wigner or Poisson curves while their av- 
erage, over the three BC, cannot fit to the semi-Poisson, 
either. However, the average over Dirichlct (hard wall) 
in both directions (00) and periodic in one direction with 
Dirichlet in the other (10), which is displayed in Fig. 3 
for various system sizes, gives a distribution very well de- 
scribed by the semi-Poisson curve of Eq. 1. This is the 



most important result of the paper and shows the va- 
lidity of the semi-Poisson for the chosen specific average 
over BC at criticality. The obtained semi-Poisson is also 
in agreement with recent results for the critical P{S) at 
the metal-insulator transition in 3Z) disordered systems, 
where the semi-Poisson was obtained by averaging over 
all possible combinations of BC pffl. 
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FIG. 3. The mean P(S) distribution of the (10) and (00) com- 
binations of BC for several system sizes L is shown to follow closely 
the semi-Poisson distribution P{S) = (3125/24)5** exp(-5S) (Eq. 
(1) for /3 = 4) (black line). The Wigner (dotted line) and the 
Poisson (dashed line) are also plotted. 

The longer, in the i?-range, critical spectral fluctua- 
tions are described by a linear number variance E2 {E) « 
Xo + xE!, with the compressibility x related to the crit- 
ical wavefunction dimension D2 and the space dimen- 



sion d via X = (1/2) (1 



Dt/d) 



model previous studies gave D2 



L6[. In the considered 
1.63 @. For a rather 
small energy window E (see Fig. 4) the level number 
variance is shown to be linear with level compressibil- 
ity X which varies with the chosen BC. However, when 
the energy window E increases the result becomes inde- 
pendent of BC, as was already shown in 3D |10|. The 
obtained asymptotic level compressibility in this case be- 
comes T,2iE)/E ^ X « 0.07 which leads to a D^, rather 
close to the expected value according to the previous for- 
mula. The obtained value of x in 2D should be contrasted 
with the higher 3D asymptotic value x ~ 0.27 pc| ]. 

The main result from our calculations, done on a 
square random lattice and not on a peculiarly shaped 
non-random billiard, is the validity of the semi-Poisson 
statistics at the metal-insulator transition in 2D disor- 
dered systems with spin-orbit coupling. However, the 
semi-Poisson P{S) is obtained for the averaged distribu- 
tion over two specific BC. It must be pointed out that 
only the main part of the distribution agrees surprisingly 
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well with the analytical result. For large S, where the de- 
pendence on BC should become less important, we have 
not suceeded to describe its tails, possibly due to their 
exponentially small nature. In this case the appropriate 
statistical measure becomes the number variance since 
longer range level correlations are needed. We find a 
linear number variance 'S2{E) ~ xE which becomes in- 
dependent of the BC choice. The obtained x is close to 
the expected value from the formula via D2 ■ 




FIG. 4. The number variance T:2{E) vs the energy window E 
for the two BC (00) and (10). The straight lines fit the data giv- 
ing slopes corresponding to non-asymptotic x values. The Wigner 
(dotted line) and the Poisson (dashed line) are also displayed for 
comparison. 

In conclusion, we have shown the validity of the semi- 
Poisson level statistics at the critical point of the metal- 
insulator transition with /3 = 4 in 2D. The semi-Poisson 
curve is shown to describe very well the main part of 
the computed distribution for a specific average over BC 
and is similar to recent results for critical disordered sys- 
tems and weakly chaotic quantum systems. Our cal- 
culations, on one hand, can give a justification for the 
averaging over boundary conditions recently shown to 
lead to the semi-Poisson statistics at the mobility edge 
1^ . On the other hand, suggest that such an average 
might be related to the bandwidth distributions, by re- 
peating periodically the square, as it was recently shown 
for a non-random one-dimensional critical quasi-periodic 
model ||l7[| . Clearly, more work is needed to examine the 
validity of the critical semi-Poisson distribution, also for 
3Z) disordered systems with spin-orbit coupling and sys- 
tems in the presence of a magnetic field {(3 = 2). 

This work was supported by a TMR network of EU and 
a Greek- Chineese collaboration G.S.S.T. Useful discus- 
sions with J.-L. Pichard and G. Montambaux are grate- 
fully aknowledged. 
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FIG. 5. Ti2{E)/E versus E for a much broader range of E 
where the independence on BC is seen. The asymptotic value ap- 
proximates X ss 0.07. 
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